Design of Strict Control- Lyapunov Functions 
for Quantum Systems with QND Measurements * 



H. Aminif P. Roucnon* M. Mirrahimi § 
January 14, 2013 

Abstract 

We consider discrete-time quantum systems subject to Quantum Non-Demolition 
(QND) measurements and controlled by an adjustable unitary evolution between 
two successive QND measures. In open-loop, such QND measurements provide a 
non-deterministic preparation tool exploiting the back-action of the measurement 
on the quantum state. We propose here a systematic method based on elementary 
graph theory and inversion of Laplacian matrices to construct strict control-Lyapunov 
functions. This yields an appropriate feedback law that stabilizes globally the system 
towards a chosen target state among the open-loop stable ones, and that makes in 
closed- loop this preparation deterministic. We illustrate such feedback laws through 
simulations corresponding to an experimental setup with QND photon counting. 

1 Introduction 

Feedback stabilization of quantum states is closely related to the concept of Quantum 
Non-Demolition (QND) measurement [3j [121 03] • Indeed, as soon as we are interested in 
applying a measurement-based feedback to stabilize a quantum state, we need to make 
sure that the measurement itself is not changing the desired target state. This means that 
the measurement procedure is QND with respect to the projection over the target state. 
In fact, very often a well-chosen QND measurement protocol can itself be considered as a 
preparation tool for various quantum states. However, this preparation is generally non- 
deterministic and one can not make sure to converge towards the desired state except by 
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repeating the experiment many times. The feedback can be applied here to make this 
process deterministic [TT| |8| II]. 

This paper is a generalization of the feedback law, proposed in [U E] and experimentally 
tested in j9], to generic discrete-time quantum systems where, between two successive 
QND measurements, a controlled unitary evolution can be applied. The dynamics of such 
discrete-time quantum systems are governed by non linear Markov chains. In jU [7] the 
feedback laws were obtained by maximizing the fidelity with respect to the target state at 
each time-step: this means that the feedback strategy was based on the same Lyapunov 
function given by the fidelity between the current and the target state. We propose here a 
systematic and explicit method to design a new family of control-Lyapunov functions. The 
main interest of these new Lyapunov functions relies on the crucial fact that the increases of 
their expectation values at step k + 1 knowing the state p k at step k remain strictly positive 
when p k does not coincide with the target state. In closed-loop, these Lyapunov functions 
become strict and the convergence analysis is notably simplified since invariance principles 
are not necessary. The construction of these strict Lyapunov functions is based on the 
Hamiltonian H underlying the controlled unitary evolution and relies on the connectivity 
of the graph attached to H . They are obtained by inverting a Laplacian matrix derived 
from H and the quantum states that are untouched by the QND measurements. 

In section [2j we describe the finite dimensional Markovian model together with the 
main modeling assumptions. Section [3] is devoted to the open- loop behavior (Theorem 13. ip 
that can be seen as a non-deterministic protocol for preparing a finite number of isolated 
and orthogonal quantum states. In Section HJ we present the main ideas underlying the 
construction of these strict control-Lyapunov functions W e . Then we define the connec- 
tivity graph, the Laplacian matrix attached to H and two technical lemmas used during 
this construction. Finally Theorem 14.11 describes the stabilizing feedback derived from 
W e . Closed-loop simulations corresponding to an experimental setup at Ecole Normale 
Superieure are sketched in section |5j 

The authors thank M. Brune, I. Dotsenko, S. Gleyzes, S. Haroche and J.M. Raimond for 
enlightening discussions and references. Advices of L. Praly concerning control-Lyapunov 
functions are also acknowledged. 

2 The non- linear Markov model 

We consider a finite dimensional quantum system (the underlying Hilbert space "H = C d is 
of dimension d > 0) being measured through a generalized measurement procedure placed 
discretely in time. Between two measurements, the system undergoes a unitary evolution 
depending on a scalar control input u£l. The dynamics of such discrete time quantum 
systems is described by a non-linear controlled Markov chain whose structure is derived 
from quantum physics. We just sketch here this structure with a mathematical viewpoint. 
A tutorial physical exposure can be found in [5]. 

The system state is described by the density operator p belonging to V(7-L) the set of 



positive, Hermitian matrices of trace one: 

V{H) := {p G C dxd | p = p\ Tr(p) = l, p>0}. 

The generalized measurement procedure admits m > different discrete values p G 
{1, . . . , m}: to each measurement outcome p is attached a Kraus operator described by a 
matrix M M G C dxd . The Kraus operators (M M ) Mg { lv .. im } satisfy the constraint 5^«Li M^M^ = 
1 where 1 is the identity matrix. In general the M M are not necessarily Hermitian. 
The controlled evolution between two measures U u is defined by the unitary operator 
exp(-iuH) = U u where if is a Hermitian operator H G <C dxd with W = H . 

The random evolution of the state pk G V(H) at time step k is modeled through the 
following Markov process: 

Pk+1 = V Uk (M, k (p k )), (1) 

where 

• Uk G R is the control at step k, 

• pk is a random variable taking values p in {1, • ■ ■ , m} with probability p PPk = 
Tr(M M p fc Mt), 

• U u is the super-operator 

U u : V(U)3 p^U uP UteVCH), 

• For each p, M M is the super-operator 

ML : p M. G £>(ft) 

defined for p G £>("H) such that Tr (M^pMj) ^ 0. 

We suppose throughout this paper that the two following assumptions are verified by the 
system under consideration. 

Assumption 1. The measurement operators M M are diagonal in the same orthonormal 
basis { \n) \ n G {1, • • • , d}}, therefore M p = Ylt,=i c n,n \ n ) ( n l w ^ c m,« e ^• 

Assumption 2. For all n x ^ n 2 in {1, • • • ,d}, there exists a p G {1, • • • ,m} such that 
Ic I 2 \c I 2 

Assumption [1] means that the considered measurement process achieves a Quantum 
Non Demolition (QND) measurement for the physical observables given by orthogonal 
projections over the states (|n) | n G {1, • • • , d}\. This implies that for u k = 0, any p = 
\n) (n| corresponding to the orthogonal projector on the basis vector |n), n G {1, . . . , d}, is a 
fixed point of the Markov process (pQ). Since the operators M M must satisfy X]jT=i M^Mfj, = 
1, we have, according to assumption [U ^™ =1 \ c v,n\ 2 = 1 for all n G {1, • • • , d}. 

Assumption [2] means that there exists a /x such that the statistics when u k = for 
obtaining the measurement result /i are different for the fixed points |ni) (ni| and \n2) {u2\- 
This follows by noting that Tr (M M |n) (n| M^) = |c Mjn | 2 for n G {1, • • • , d}. 



3 Convergence of the open-loop dynamics 

When the control vanishes (u k = 0, Vfc), the dynamics is simply given by 

p k+1 =M Pk (p k ) (2) 

where //& is a random variable with discrete values in {1, . . . , m}. The probability p^ Pk to 
have //fe = \i depends on p k : p PtPk = Tr (M M pfcM£). We have then the following theorem 
characterizing the open-loop asymptotic behavior: 

Theorem 3.1. Consider a Markov process p k obeying the dynamics of with an initial 
condition p in Then 

• with probability one, p k converges to one of the d states \n) (n\ with n e {1, • • • , d}. 

• the probability of convergence towards the state \n) (n| depends only on the initial 
condition po and is given by 

Tr{p |n) (n|) = (n\p \n) . 

The proof is a generalization of the one given in p] . 
Proof. For any n G {1, . . . , d}, Tr (|n) (n| p) is a martingale. This results from 

E(Tr(|n) (n\ p k +i) \pu) 

m 

= ^ Tr (M M p fc Mt) Tr (|n) (n| M„(p fc )) 



/x=l 

/ m 

^ " ' r 1 'n| p fc 



£ Tr (|n> (n| M M p fc Mt) = Tr PT |n) ( 



Tr(|n) (n|p fc ) 



where we have used the facts that M p and |n) (n| commute and that Y2™=i M^M^ = 1. 
Take the function 

d 

V(p):=^/C&(|n)<n|p)). (3) 

n=l 

where f(x) = -|-. The function / being convex and each Tr (|n) (n| p) being a martingale, 
we infer that V(p) is a sub-martingale 

E(y(p fc+ i)|p Jfc ) > v( Pk ). 

More precisely, we have 
E(V(j> k+1 )\p k ) = 

EE 11 (^m.) / 



n=l /tG/, 



with I Pk = jp G {1, • " " > m } I Tr (M^p^M^) ^ oj. For all sequence of reals Xi, • • • ,x m 
and $i, ■ ■ ■ , 9 m in the interval [0, 1] with Y^=i ®n = 1> we have the identity 

m / m \ m 

E e ^ w = / E ^ + E w***^' ( 4 ) 

For each p, let 6^ = Tr (M^p&Mj), a; M = ~^7^~^F^ if 6^ > and x M = otherwise. 
Identity © yields 



E(V(p fc+ i)|pfc) - V(pfc) 



■1 



lE E Tr (M M p fe Mt) Tr (M,p fc Mt) . . . 
n=l n,uei Pk 

\ 2 

|e M , n | 2 (n|pfc|n) |c„, n | 2 (n|p fc |n) 



(5) 



We have used the fact that (n| M M (p fe ) |n) = ^t^tj - Thus, V(p fc ) is a sub-martingale, 

and in addition we have a precise bound on the difference E (V(pk+i)\pk) — V(Pk)- We 
denote by Q{pk) the right term in Equation (jSj). Note that the sum in the definition of 
Q{pk) is over all p, z/ G J Pfc . However, we can assume that the sum is actually over all p, z/ 
in {1, • • • , m} by observing the following facts. 
For any p, u, take the mapping 



Tr f M oM ] \ Tr f M oM^ I l c ^"l 2 ( re l^) _ Kn| 2 Hp|n) 



defined only when Tr (M M pMj) Tr (M u pMl) > 0. Since this mapping is positive and 
bounded by p Tr (M^pMj) Tr (M u pM^, it can be extended by continuity to any p G 
by taking a null value when Tr (M^pMj) Tr (M v pMl) = 0. Thus 



(6) 



w = i E E Tr ( M ^K) Tr i M »p M D 

n=l p,z/=l 

\ 2 

\cn,n\ 2 {n\p\n) _ |c„ iT1 | 2 (n|p|n) 
Tr(M M pA/t) Tr(A^pMj) 

is continuously defined for any p G T>(H) and still satisfies 

E(y(p fc+1 )|p fe ) - V( Pk ) = Q( P k),V Pk G X?(H). 

By Theorem IA.1I of the Appendix, the w-limit set Q (in the sense of almost sure conver- 
gence), for the trajectories pk, is a subset of the set {p G V{%)\ Q(p) = 0}. 



Let us consider a density matrix p^ in the w-limit set. Therefore Q(poo) = implies 

for all p and v in J Poo and for each basis element \n). Since Tr (poo) = 1, there is at least 
one n such that (n| p^ |n) > 0. Then Equation ([7]) simplifies to 

Tr (M^ Poo Ml) \c u , n \ 2 = Tr (M vPoo Mt) |c^| 2 , 

for all p, z/ G J Poo . Summing over z/ G J Poo , we find 

2 X 



Tr {M, Poo Ml) ( ]T |c^| 



\c -\ 2 ) 






1 C[i,n | 



( ^ Tr (M vPoo Mt)) | 

By definition of I Poo 

m 

Tr (M vPoo Ml) = J2 Tr (M^M*) = 1. 

Thus we have 

Tr (M MPoo Mt) ( 

For z/ J Poo , = since X)n=i l c ^,«| 2 ( n l Poo \n) = 0, each (n| Poc \n) > and (n| p^ |n) > 
0. Thus S^/^ l c ^| 2 = XT=i \ c vM 2 = !■ Finally, we have 

Vp G I Poo , Tr (M pPoo Mt) = |c^| 2 . (8) 

as soon as (n \ Poo \n) > 0. 

Assume now that exist n\ ^ n-i in {1, . . . , d) such that (n\ \ Poo \ni) > and (rz 2 1 Poo 1^2 ) > 
0. Then (jSJ) implies that 

Vp G / poo , | C At,ni| = | c p,n2l • (9) 

By Assumption [21 there exists a p G {1, ■ • • , m} such that |c p>Sl | 2 ^ |c Pi n 2 | 2 - These terms 

cannot be simultaneously zero, thus Tr yM^p^Mp^j > 0, p G I Poo . This is in contradiction 

with (J9j). This closes the proof of the assertion: the w-limit set is reduced to the set fixed 
point \n) (n\ , with n G {1, • • • , d}. 

We have shown that the probability measure associated to the random variable pk 
converges to the probability measure ^ n= iPn^|n>(n|> where 6\ n )( n \ denotes the Dirac mea- 
sure at |n) (n\ and p n is the probability of convergence towards \n) (n\ . In particular, we 
have E (Tr (p*. \n) (n\)) — > p n . But Tr (|n) (n\ pf.) is a martingale thus E (Tr (|n) (n\ p^)) = 
E (Tr (\n) (n| po)) and consequently p n = (n\ p \n). □ 



4 Feedback stabilization 



4.1 Design of strict control-Lyapunov functions 

The goal is to design a feedback law that globally stabilizes the Markov chain ([Q) towards a 
chosen target state \n) (n\, for some n among {1, • • • ,d}. In previous publications P [TJ H], 
the proposed feedback schemes tend to increase at each step the same open-loop martingale 
Vn(p) — (n\ p\n): u is chosen in order to increase u >-)■ V s (U u (p)). When p = \n) (n\ with 
n 7^ n, u (->■ V n (p) is minimum at u — 0. Consequently, its first-order ^-derivative vanishes 
at u = 0. But its second-order ^-derivative could also vanishe at u = 0. For the photon-box 
considered in [1] (see also section [5]), this happens when |n — n| > 2. Such lack of strong 
convexity in u when the image of p is almost orthogonal to \n), explains the fact that, in 
previous works, the control u is set to a constant non-zero value when V n {p) is close to 0, 

To improve convergence and avoid such constant feedback zone, we propose to modify 
Vn using the other open-loop martingales V n (p) = (n\ p \n) and the sub-martingale V(p) = 
Ylt,=i(( n \ P\ n )) 2 use d during the proof of theorem I3.ll The goal of such modification is 
to get control Lyapunov functions still admitting a unique global maximum at \n) (n\ but 
being strongly convex versus u around when p is close to any \n) (n\, ii / n. 

Take W (p) = Y^t=i a n ( n \ P \ n ) with real coefficients a Tl to be chosen such that a n 
remains the largest one and such that, for any n ^ n, the second-order w-derivative of 
W^Uudn) (n\)) at u = is strictly positive. This yields to a set of linear equations (see 
lemma fl~2i) in a n that can be solved by inverting a Laplacian matrix (see lemma l4~Ti) . 
Notice that Wq is an open-loop martingale. To obtain a sub-martingale we consider (see 
theorem 14. lj) W e (p) = Wo(p) + eV(p). For e > small enough: W e (p) still admits a 
unique global maximum at |n) (n|; for u close to 0, W e (U u (|n) (n\)) is strongly convex for 
any n ^ n and strongly concave for n = n. This implies that W e is a control-Lyapunov 
function with arbitrary small control (see proof of theorem 14. ip . 

Let us continue by some definitions and lemmas that underlay the construction of these 
strict control-Lyapunov functions W £ . 

4.2 Connectivity graph and Laplacian matrix 

To the Hamiltonian operator H defining the controlled unitary evolution U u = e~ mH , we 
associate its undirected connectivity graph denote by G . This graph admits d vertices 
labeled by n 6 {1, . . . , d}. Two different vertices n\ ^ ri2 (ni, G {1, . . . , d}) are linked 
by an edge, if and only if, (m| H \n^) ^ 0. Attached to H, we also associate R , the real 
symmetric matrix d x d (Laplacian matrix) with entries 



Lemma 4.1. Assume the graph G H to be connected. Then for any positive reals X n , 
n G {1, . . . ,d}, n =fi n, there exists a vector a = (cr n )ne{i,-,d} °f^- d suc h that R H a = —A 
where X is the vector ofM. d of components X n for n ^ n and Xn = — Xln^n ^n- 




(10) 



Proof. Note that R is symmetric and the sum of the entries for any column and any row 
of R H is equal to zero. Therefore, the vector (1 • • • 1) T is in the kernel of R . The diagonal 
(resp. non-diagonal) components of R are positive (resp. negative). Therefore R is a 
Laplacian matrix (see [21 Ch. 4]). The connectivity graph associated to R H coincides 
with G H . Since this graph is supposed connected, classical results of graph theory (see, 
e.g., [2j Theorem 3.1]) imply that the "constant" vector (1, • • ■ , 1) T spans the kernel of R H . 
Therefore, the dimension of the image of R H is equal to d — 1. Since R H is symmetric, its 
image coincides with the orthogonal to its kernel. For the sake of completeness, here we 
give a simple proof of this statement. Indeed, for any vector X in the kernel of R , we 
have 



R 



H 

m,n 2 



X T RX 



0. 



^1,"2S{1,"- ,d} 

This implies R^ n2 (X ni — X n2 ) 2 = 0, Vrii, n 2 G {1, • • • , d}. As the graph of R H is connected, 
we necessarily have X ni = X n2 for all n 1; n 2 G {1, • • ■ , d}. Thus any vector orthogonal to 
(1 ■ ■ • 1) T is in the image of R H . The vector A is orthogonal to (1, • ■ ■ , 1) T . □ 



Lemma 4.2. Take n 6 {1, . . . , d} and consider X n > 0, for n G {1 

G H connected and consider the vector a = (cr n ) G 
p G ViTi) we set 



,d}/{n}. Assume 



W (p) = J2^Tr(\n) (n\p) 



n=l 



Then for any n G {1, . . . , d}/{n} we have 

d 2 W (v u (\n){n\)) 



u=0 



and 



d 2 W (u„(|n>(n|)) 



du 2 



u=0 



d given by Lemma \4-l\ For any 

d 

■- ^cr„ (n\ p \n) . (11) 

n=l 

Xn > 

- J2 A « < 

- W (e- mH \n) (n\e luH ). 



Proof. For any n, set g n (u) = Wo(U u (\n) (n\) = W (e %uH \n) (n\ e luIi ). The Baker- 
Campbell-Hausdorff formula yields up to third order terms in « ([•,•] is the commutator): 

U„(|n) (n|) w \n) (n| - iu[H, \n) (n|] - f [if, [If, |n) (n\}}. 
Consequently for any / G {1, . . . , d}, we have 

Tr (U u (|n) (n|) |/) (/|) » Tr (|Z) (Z| ( |n) <n| - w[IT, |n) (n|])) 

-^Tr (\l)(l\([H, [H,\n) (n\}])) 

%, n + £Tr([H,\n) (n\][H,\l) 



since Tr(|Z) (l\ [H, \n) (n 
|n) (n|), and since 



(12) 

— Tr([|Z) (l\ , |n) (n\]H) = because |Z) (l\ commutes with 
Tr (|0 (l\ [H, [H, \ n ) (n\}]) = -Tr ([H, \n) (n\][H, \l) . 



Thus up to third order terms in u, we have 



9n(u) 

Therefore: 

d 2 W (v u (\n)(n\)) 



J2<ri(5 l>n + ^Tr([H,\n) (n\][H,\l) 



1=1 



du 2 



u=0 



5>Tr([tf» (n\][H,\l) </|]) 



It is not difficult to see that Tr([if, \n) (n\][H, \l) 



1=1 



-<i Thus 



d 2 Wo(v u (\n){n\)) 



du 2 



u=0 



4.3 The global stabilizing feedback 

The main result is expressed through the following theorem. 

Theorem 4.1. Consider the controlled Markov chain of state p k obeying ([1]). Assume that 
the graph G H associated to the Hamiltonian H is connected and that the Kraus operators 
satisfy assumptions^ and\^ Take n G {1, . . . , d} and d — 1 strictly positive real numbers 



\ n > 0, n G {1, . . . , d}/{n}. Consider the component (a n ) ofa£M. d defined by Lemma 4-1 
^ii k {.Pk) the quantum state just after the measurement outcome pk at 



Denote by p k+ i - 
step k. Take u > and consider the following feedback law 



Uk = K(p k+ i) = argmax(w e (lS u (p k+ i))) : 



iiS[-u,ii] 

where the control- Lyapunov function W e (p) is defined by 

d 



(13) 



(14) 



n=l 



with the parameter e > not too large to ensure that 

Vn G {l,...,d}/{n}, X n + e (((n| H \n)) 2 - (n\ H 2 \n}) > 0. 

Then, for any po £ 'Dfti), the closed-loop trajectory p k converges almost surely to the pure 
state \n) (n\. 

The proof relies on the fact that W e (p) is a strict Lyapunov function for the closed- loop 
system. 



Proof. We have 

E(W e {p k+l )\p k )-W e ( Pk ) = 

E max (w e (U„(M M (p fc )))) - W e (p fc ) 

E ^, Pfc (^(M /t (p fc ))-^(p fc )) + 



AteJ 



We define the following functions of p k 



E P», Pk ( max (W £ (U U (M^(^)))) - W e (M M (p fc ))). 



and 



E max (w e (U„(M> fc )))) -W e (M> & ))). 



These functions are both positive continuous functions of p k (the continuity of these func- 
tions can be proved in the same way as the proof of the continuity of Q{pk) in Theorem l3.ip . 
By Theorem IA.1I of the appendix, the w-limit set Q is included in the following set 

{ P e V{U)\ Qi( P ) = 0} n { P e v(U)\ Q 2 (p) = o}. 

Indeed Qi coincides with Q defined in ([6]). During the proof of Theorem I3.1[ we have 
shown that Qi(p) = implies p = \n) (n\ for some n G {1, • • • , d} . But Q%{\n) (n\) = 
implies that max W e (U u (\n) (n\)) = W e ({n\ \n) ) since M M (|n) (n\) = \n) (n\. According 



to the Lemma [4.21 and the relation ( TT2l) . we have We ^ V "jl ~ 



= and 

u=0 



d 2 W € (v u (\n){n\)) 



du 2 



= X n + e({(n\H\n)) 2 - (n\H 2 \n)) 

u=0 



Since A n > for n ^ n and e is not too large, for any n ^ n, u = is a locally strict 
minimum of W e (U u (|n) (n\)). Consequently Qzdn) (n|) = implies that n = n. □ 



5 Closed-loop simulations for the Photon-Box 



We recall the photon-box model presented in p£]: d = n max + 1, where n max is the maximum 
photon number. To be compatible with usual quantum-optics notations the ortho-normal 
basis |n) of the Hilbert space H = C" max+1 is indexed by n e {0, . . . , n max } (photon number). 
For such system p k takes just two values g or e, and the measurement operators M g and 
M e are defined by M g = cos(0o + #N) and M e = sin(0 o + #N) (</>q and 6 are constant 
angles fixed as the experiment parameters). When 9 /it is irrational, assumption [5] is 
satisfied. The photon number operator N is defined by N = a^a, where a is the annihilation 
operator truncated to n max photons, a corresponds to the upper 1-diagonal matrix filled 
with (vT, ■ ■ • , \/n max ) : and N to the diagonal operator filled with (0, 1, ... , n max ). The 
truncated creation operator denoted by at is the Hermitian conjugate of a. Consequently, 
the Kraus operators M g and M e are also diagonal matrices with cosines and sines on the 
diagonal. 

The Hamiltonian H = z(a^ — a) yields the unitary operator U u = e u(at ~ a ) (also known 
as displacement operator). Its graph G H is connected. The associated Laplacian matrix 
R H admits a simple tri-diagonal structure with diagonal elements R^ n = 4n + 2, upper 
diagonal elements R^-i n = — 2n and under diagonal elements R^+x n = ~ 2 (up to 
some truncation distortion for n = n max — l,n max ). 

For a goal photon number n G {0, . . . ,n max — 1}, we propose the following setting for 
the X n and e defining W e of Theorem 14.11 

e {0,...,n max }/{n}, A n = 1. 

Figure [1] corresponds to the values of o n found with X n given in above with n = 3 and 
n max = 10. We remark that a n is maximal. 

Since ((n\ H \n)) 2 — (n\ H 2 \n) = —2n — 1, the constraint on e > imposed by Theo- 
rem O reads Vn ^ n, e < i.e., e < 2nm L +r 

The maximization defining the feedback in Theorem 14.11 could be problematic in prac- 
tice. The unitary propagator U u does not admit in general a simple analytic form, as it is 
the case here. Simulations below show that we can replace U u by ist quadratic approxima- 
tion valid for u small and derived from Baker-Campbell-Hausdorff: 

U u (p) =p- iu[H, p] - f[H, [H, p}) + 0(\u\ s ). 

This yields to an explicit quadratic approximation of W e (JJ u (p)) around 0. The feedback 
is then given by replacing W e (V u (p)) by a parabolic expression in u for the maximization 
providing the feedback Uk = K(p k+ i). 

For the simulations below, we take n max = 10, n = 3 and e = , m L , an d u = 

_ 4n max +2 10 

The parameters appearing in the Kraus operators are 9 = ^ and 0o = f — n9. Figure |2] 
corresponds to 100 realizations of the closed loop Markov process with an approximated 
feedback obtained by the quadratic approximation of U u versus u sketched here above. 
Each realization starts with the same initial state U /^(|0) (0|) (coherent state with an 
average of 3 photon). Despite the quadratic approximation used to derived the feedback, 
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Figure 1: The coefficients a n used for the Lyapunov function W e (p) defined in Theorem 14.11 
during the closed-loop simulations, the goal photon number n = 3 corresponds to the 
maximum of the a n . 

we observe a rapid convergence towards the goal state |n) (n|. In more realistic simulations 
not presented here but including the quantum filter to estimate pk from the measurements 
Pk and also the main experimental imperfections described [I] , the asymptotic value of the 
average fidelity is larger than 0.6. This indicates that such feedback laws are robust. 

6 Concluding remarks 

The method proposed here to derive strict control-Lyapunov could certainly be extended 
to 

• prove exponential closed-loop convergence for the feedback law given by theorem 14. 1[ 

• the general situation where the control u appears directly in the Kraus operators 
M^{u) instead of being separated from the QND measures and attached to a unitary 
evolution applied after each measurement. 

• continuous-time quantum systems subject to QND measurements such as those con- 
sidered in [8] and described in detail in [T4"] . 

• infinite-dimensional quantum systems as in [TO] that consider the photon-box system 
without truncation to a finite number of photons. 
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A Appendix 

Theorem A.l. Let X k be a Markov chain on the compact state space S. Suppose, there 
exists a non-negative function V(X) satisfying 

E(V(X k+1 )\X k ) - V(X k ) = Q(X k ), (15) 

where Q{X) is a positif continuous function of X, then the cu-limit set Q (in the sense of 
almost sure convergence) of X k is included in the following set 

I:={X\ Q(X) = 0}. 



Proof. The proof is just an application of the Theorem 1 in [SI Ch. 8], which shows that 
Q(X k ) converges to zero for almost all paths. It is clear that the continuity of Q(X) with 
respect to X and the compactness of S implies that the w-limit set of X k is necessarily 
included into the set I. □ 



